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Introduction 


Several approaches have been used to estimate the elastic 
lithosphere thickness (L) on Venus. In particular, Solomon and 
Head C13 have used a number of compressional (elastic and 
viscous folding) and tensional (wedge subsidence, imbricate 
normal faulting, and plastic necking) models to predict L based 
on the hypothesis that the linear bands of greater and lesser 
backscatter observed in earth-based radar images of Ishtar Terra 
are of tectonic origin. All of these models predict very thin 
lithospheres, on the order of 0.3 to 8 km. These values are in 
agreement with those obtained from strength calculations based 
on laboratory measurements of crustal rocks (C13; see below) 
which give negligible yield stresses at shallow depths due to 
the high surface temperature and a probable earth-like thermal 
gradient C23. 

However, rift-like features which are at least 
superficially similar in size, morphology, and association with 
volcanism and doming to continental rifts on the earth C3D also 
exist on Venus. They typically have widths on the order of 
75-100 km. Estimates of L derived from the widths of these 
features in Aphrodite Terra and Beta Regio using a 
wedge-subsidence model are on the order of 50-70 km C 43, which 
is inconsistent with our knowledge of rock rheology at probable 
Venus temperatures. Plastic necking models also require a 
relatively thick, strong high viscosity zone (about 30 km thick 
C53), again too thick for expected crustal temperatures. 

How can these large rift structures which appear to require 
a thick, strong lithosphere be reconciled with the hot, 
apparently thin lithosphere implied by rock rheology at venusian 
conditions and the spacing of the banded terrain? In this 
abstract we briefly review lithospheric strength envelopes and 
explore their implications for large scale rifting on Venus. We 
then use these results to constrain possible crustal thicknesses 
and thermal gradients. 

Lithospheric Strength Envelopes 

The maximum stress levels found in the earth's crust are 
accurately predicted by Byerlee's law C6,73. This relation is 
based on laboratory measurements of the frictional resistance to 
sliding on pre-existing fractures, which occurs at stresses less 
than those required to break intact rock. Byerlee's law is of 
the form a^yo^+b, where a and a are, respectively, the 
maximum and minimum principle effective 3 stresses (stress minus 
pore pressure), y is the coefficient of friction, and b is a 
constant. Laboratory friction measurements show that y and b 
are virtually independent of stress (except for a slight change 
at 135 MPa), rock type, displacement, surface conditions, and 
temperature. Because the vertical stress is generally quite 
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close to the lithostatic load, this relation predicts a linear 
increase in yield stress with depth. 

With increasing temperature, rock deformation occurs 
predominantly by ductile flow. Flow laws for many rocks and 
minerals have been experimentally_^etermined for stresses up to 
1-2 GPa and strain rates down to 10 /sec. These results can be 
extrapolated to geological strain rates n via creep equations, 
which generally are of the form s=A(a -a )exp( -Q/RT) , where e 
is the strain rate, R is the gai Constant, T is absolute 
temperature, and A, Q (the activation energy), and n are 
experimentally determined constants. As a result, the ductile 
strength is negligible at depths where T is high and increases 
exponentially with decreasing depth. Flow laws are also highly 
dependent on rock composition, with a silicic crust much weaker 
than a mafic mantle. 

The failure criterion for a given depth in the lithosphere 
is determined by the weaker of the frictional or ductile 
strength at that depth (Figure 1). The yield stress increases 
with depth according to Byer lee's law until it intersects the 
crustal flow law. It then decreases exponentially until it 
reaches the moho, where the mantle flow law causes an abrupt 
increase in yield stress followed by another exponential 
decrease. In actuality, semibrittle and low temperature ductile 
processes tend to round off the intersection points between the 
brittle and ductile curves in Figure 1 C83. The lithosphere 
will behave elastically for stresses less than the yield 
stress . 

This type of analysis has been used by many investigators 
to study lithospheric strength on the earth and the terrestrial 
planets (e.g. L9-123). However determining the geometry and 
characteristics of faults at failure for this type of 
lithospheric model has not received much attention. Based on an 
analysis of two phases of extension in the Rio Grande rift, 
Morgan and Golombek L133 suggested that during the early phase a 
shallow brittle-ductile transition along with the absence of a 
zone of upper mantle strength (due to an elevated geotherm) 
allowed a shallow decollement to develop with significant strain 
between the brittle upper crust and the ductilely extending 
material below. Thus, large strains resulted in the strongly 
rotated blocks bounded by numerous listric (curved) or planar 
normal faults which are characteristic of this period (see also 
Smith and Bruhn C103). During the late phase of extension that 
formed the present horst and graben physiography the geotherm 
had cooled enough to allow a significant zone of ductile 
strength in the upper mantle (although most of the strength was 
still in the crust) which prevented large scale intracrustal 
decoupling and steep rotation of upper crustal blocks. 

Lithospheric Strength on Venus 

Figure 1 illustrates a typical lithospheric strength curve 
for Venus with a crustal thickness (c) of 10 km. For these 
calculations we assume a crustal density (P c ) of 2.8 Mg/m , a 
mantle density (p ) of 3.3 Mg/m 3 , a surface temperature of 
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740°C, a thermal gradient of 12 /km C23, and e=10 /sec 
(equivalent to about 3% extension per million years). Based on 
surface geochemical measurements which indicate a basaltic 
composition, we use a dry diabase flow law C143 for the crustal 
layer. A dry olivine flow law C153 is assumed for the mantle. 
In this abstract we are interested in horizontal tensional 
stresses, which correspond to the left-hand branch of the curves 
in Figure 1. 


Implications for Rift Width and Crustal Thickness 

One consequence of this model is that the mechanical 
structure of the lithosphere depends strongly on the crustal 
thickness. For crusts thicker than about 30 km, there is no 
appreciable strength in the mantle and the brittle extension of 
the upper crustal layer should be effectively decoupled from the 
mantle. This was the implicit assumption of Solomon and Head 
C13. However if the crust is thin, the lithospheric strength 
will be dominated by the upper mantle, and it is possible that 
this layer will control the style of deformation. In this case 
the effective lithosphere which controls rift width would be 
much thicker than the few km expected from crustal rock 
rheology. Instead, the elastic lithosphere thickness would be 
on the order of 20-30 km, based on an olivine flow law. 

Here we will consider the implications for one theory for 
rift formation, the wedge -subsidence hypothesis of Vening 
Meinesz C16,173. This theory predicts that the graben, width w 
will be given by w=ira./4, where a=(EL /3gp(l-v )) , E is 
Young's modulus, g is gravitational acceleration, p is the 
difference in density between the layers above and below the 
faulted layer, and v is Poisson's ratio. For rifting of a 
crustal layer, p is just p . However for a layer at depth with 
an effectively fluid laye? above it, p=p -p . Assuming 
E=l. 25x10 Pa, v= . 25 , p -p =0.5 Mg/m , and E=20 km gives the 
width for a graben formed iP? tfte mantle brittle zone of 98 km, 
in good agreement with observed rift widths on Venus' surface. 

This result requires that the crust be no thinner than ^3 
km (in order for there to be a viscous layer over the mantle 
strong zone) and no thicker than ^20 km (in order for there to 
be a brittle zone in the mantle). In addition, the thermal 
gradient could not be significantly greater or less than 
10-15° /km, since L should be ^20 km. It should be noted that 
these values depend on flow laws which are somewhat uncertain 
due to their extrapolation from laboratory conditions. 

This model is not inconsistent with the thin lithosphere 
implied by the banded terrain. One might expect smaller scale 
fault blocks whose widths are controlled by the thickness of the 
crustal strong layer to be superimposed upon the rift structure 
(see Figure 2). In Ishtar Terra, where no rifts are observed, 
the bands could have been formed by compressive stresses Cl, 183, 
or the crust in this region might be thicker than in rifted 
areas C2,183. 
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Figure 1 . Lithospheric strength envelope for Venus. Stress 
difference is the vertical minus the maximum or minimum 
horizontal stress. Thus negative values denote tension and 
positive values correspond to compression. The horizontal 
line at 10 km shows the crust-mantle boundary assumed for 
this case. The flow laws used for the crust and mantle are 
dry diabase C14D and oleine C151, respectively; both 
assume a strain rate of 10 /sec, a Q surface temperature of 
730 C, and a thermal gradient of 12 /km. The properties of 
the layers delineated by the tensional envelope are shown 
to the left. The elastic /ductile portion of the mantle 
deforms ductilely when its yield stress is exceeded, but 
behaves elastically at the stress levels required for 
ductile flow in the lower crust. 





Figure 2 . Schematic diagram of a modified Vening Meinesz model 
for rift formation on Venus, with the strong mantle layer 
controlling rift width and the brittle crustal layer 
controlling fault spacing. 
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